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Abstract 

The first measurement of inclusive Djr photoproduction at HERA has been per- 
formed with the ZEUS detector for photon-proton centre-of-mass energies 130 < 
W < 280 GeV. The measured cross section for 3 < p^ s < 12 GeV and \r] Ds \ < 1.5 is 
(Jep^D a x = 3.79 ± 0.59 (stat.) 1^46 (syst.) ± 0.94 (br.) nb, where the last error arises 
from the uncertainty in the Djr decay branching ratio. The measurements are 
compared with inclusive D*^ photoproduction cross sections in the same kinematic 
region and with QCD calculations. The D^ 1 cross sections lie above a fixed-order 
next-to-leading order calculation and agree better with a tree-level 0(ac$.) calcu- 
lation that was tuned to describe the ZEUS D* ± cross sections. The ratio of 
to D*± cross sections is 0.41 ± 0.07 (stat.) tom (syst.) ± 0.10 (br.). From this ratio, 
the strangeness-suppression factor in charm photoproduction, within the LUND 
string fragmentation model, has been calculated to be 7 S = 0.27 ± 0.05 ± 0.07 (br.). 
The cross-section ratio and 7 S are in good agreement with those obtained in charm 
production in e + e~ annihilation. 
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1 Introduction 



Inclusive Df photoproduction cross sections at HERA are presented for photon-proton 
centre-of-mass energies in the range 130 < W < 280 GeV. The Df mesons were recon- 
structed through the decay chain Df — ► 07r ± — ► (K + K~)n :iz . This analysis supplements 
recent measurements of inclusive photoproduction of D*^ mesons at HERA |l], [J. The 
high-statistics measurement by the ZEUS collaboration [TLJ was performed in the same 
W range as given above. The measured cross sections were compared to next-to-leading 
order (NLO) calculations [|, [5| with fragmentation parameters extracted from D*^ 
production in e + e~ annihilation. The experimental results were found generally to lie 
above the NLO expectations, in particular in the forward (proton) direction. 

The study of Df photoproduction provides another test of perturbative QCD (pQCD) cal- 
culations of charm production |3], [(J which is experimentally independent of the D*^ mea- 
surement. Furthermore, from a ratio of the Df and D*^ cross sections the strangeness- 
suppression factor, 7 S , in charm fragmentation can be determined. A comparison of the 
cross-section ratio and 7 S with those obtained in charm production in e + e~ annihilation 
tests the universality of charm fragmentation. 



2 Experimental Conditions 

The measurements were performed at the HERA ep collider in the ZEUS detector during 
1996/1997. In this period HERA collided positrons with energy E e = 27.5 GeV and 
protons with energy E p = 820 GeV. The integrated luminosity used in this analysis is 
38 pb^ 1 . A detailed description of the detector can be found elsewhere f7j. 

Charged particles were measured in the central tracking detector, CTD ||, which is 
a drift chamber consisting of 72 concentric sense-wire layers covering the polar angle?] 
region 15° < 9 < 164°. The CTD operates in a magnetic field of 1.43 T provided by a thin 
superconducting solenoid. The transverse momentum resolution for full length tracks is 
v P± /p± = 0.0058 p ± ® 0.0065 ©0.0014/p_j_ (p ± in GeV). To estimate the energy loss, dE/dx, 
of tracks, the truncated mean of the sense-wire pulse-heights was recorded for each track, 
discarding the 10 % lowest and 30 % highest pulses. 

The solenoid is surrounded by the uranium-scintillator sampling calorimeter (CAL) ||, 
which is almost hermetic and consists of 5918 cells, each read out by two photomultipliers. 
Under test beam conditions, the CAL has a relative energy resolution of 0.18/v^E (E in 
GeV) for electrons and 0.35 /y/E for hadrons. 

The luminosity was determined from the rate of the bremsstrahlung process e + p — > e + 7p, 
where the photon was measured by a lead-scintillator calorimeter [ID located at Z — 
-107 m. 



1 The ZEUS coordinate system is right-handed and has the nominal interaction point at X = Y = Z = 
0, with the Z-axis pointing in the proton beam direction and the X-axis horizontal, pointing towards the 
center of HERA. 
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3 Event Selection and Reconstruction 

s 

The ZEUS detector uses a three-level trigger system [[/J. At the first level (FLT), the 
calorimeter cells were combined to define regional and global sums that were required to 
exceed various CAL energy thresholds. At the second level (SLT), beam-gas events were 
rejected by cutting on the quantity Hi(E — pz)% > 8 GeV, where the sum runs over all 
calorimeter cells and pz is the Z component of the momentum vector assigned to each 
cell of energy E. At the third level (TLT), at least one combination of CTD tracks was 
required to be within wide mass windows around the nominal (f) and D s mesonP] mass 
values, assuming ir and K masses as appropriate. The D s transverse momentum, p® s , 
was required to be greater than 2.8 GeV. 

Photoproduction events were selected by requiring that no scattered positron was identi- 
fied in the CAL ]TJ. The Jacquet-Blondel estimator of W, W m = y/2E p Xi(E - p z )i, 



was required to be between 115 and 250 GeV. The lower limit was due to the SLT re- 
quirement, while the upper one suppressed remaining DIS events with an unidentified 
scattered positron in the CAL flllf . After correcting for detector effects, the most im- 
portant of which were energy losses in inactive material in front of the CAL and particle 



losses in the beam pipe |TT| , [13|| , this Wjb range corresponds to an interval of true W 
of 130 < W < 280 GeV. Under these conditions, the photon virtuality, Q 2 , is limited 
to values less than lGeV 2 . The corresponding median Q 2 was estimated from a Monte 



Carlo (MC) simulation to be about 3 x 10 GeV . 

The MC sample used for this analysis was prepared with the PYTHIA 6.1 |14|| gener- 



ator. The proportions of direct- and resolved-photon events JOJ corresponded to the 
PYTHIA cross sections (~ 50% each), where charm excitation processes were included 
in the resolved component M. The MRSG [[LlJ and GRV-G HO [ff7 parametrisations 



were used for the proton and photon structure functions, respectively. The MC events 
were processed through the standard ZEUS detector- and trigger-simulation programs 
and through the event reconstruction package used for offline data processing. The data 
and MC distributions were found to be in good agreement. 



The D s mesons were reconstructed through the decay mode Df — > (pir^, which has a 



branching ratio Bp^^ = 0.036 ± 0.009 The was identified via its decay mode 



(f) -> K+K~, with 

^k+k- — 0.491 ± 0.008 ||18|| . The analysis was restricted to the 
pseudorapidity range —1.5 < r] Ds < 1.5, for which the CTD acceptance is high. Here 
r] Ds = — ln(tan |), where 9 is the polar angle with respect to the proton beam direction. 
The kinematic region in p^ s was limited to 3 < pf s < 12 GeV. The lower cut was required 
to comply with the p®° cut applied at the TLT. The upper cut was due to the limited 
statistics. 

The three D s decay tracks were required to originate from the event vertex, which was 
measured with a resolution of 0.4 cm in the Z direction and 0.1 cm in the XY plane. 
Only tracks with polar angles 20° < 9 < 160° and transverse momenta p± > 0.75 GeV 
were considered in the analysis. 

The decay of the pseudoscalar D s meson to the <fi (vector) plus 7r (pseudoscalar) final 
2 Here D* and D s refer to D*^ and Df, respectively. 
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state results in an alignment of the spin of the meson with respect to the direction of 
motion of the relative to the D s . Consequently, the distribution of cos9* Kl where 9* K is 
the angle between one of the kaons and the pion in the rest frame, followed a cos 2 9* K 
shape, implying a flat distribution for cos 3 9* K . In contrast, the cos6* K distribution of the 
combinatorial background was flat and its cos 3 9* K distribution peaked at zero. A cut 
of | cos 3 9* K | > 0.15 suppressed the background by a factor of approximately two while 
reducing the signal by 15%. 

The dE/dx information was used to allow partial K and 7r separation. For each kaon 
(pion) candidate a likelihood function £k(-r) = exp{ — \[{dE/dx) meas — (dE / dx) k(tt)} 2 / °' 2 } 
was determined, where (dE/dx)K(n) is the expected value for a kaon (pion), and a is 
the dE/dx resolution, which is inversely proportional to \fn, where n is the number of 
hits entering the truncated mean. The parametrisation for the dE/dx expectation was 
obtained [IS] from a fit to an independent inclusive track sample (Fig. |l]a). A normalised 



likelihood function Li = £{/ ^ £j was defined, where the sum extends over the considered 
particle hypotheses it, K and p. Provided that the number of hits was sufficiently high 
in > 7), low-likelihood hypotheses were rejected if £i < 0.05, unless Li > 0.12. The dE/dx 
cuts reduced the combinatorial background by approximately 20%. The signal loss was 
determined by means of a Monte Carlo simulation, using dE/dx parameters obtained 
from the data. The overall loss was 2.1%. 

For the D s candidates, / E e > w ° > 0.18 was required, where i?^ >10 ° is the transverse 
energy outside a cone of 9 = 10° defined with respect to the proton direction. This cut 
removed more than 20% of the background while preserving about 95% of the D s signal, 
as verified by MC studies. 

The K + K~tc ± mass distribution with the above cuts is shown in Fig. |]b for events in 
the mass range, 1.0115 < M(K + K~) < 1.0275 GeV. The fraction of events with more 
than one entry in the D s mass region was less than 1%. The mass distribution was 
fitted to a sum of a Gaussian with the D s mass and width as free parameters, and an 
exponential function describing the non-resonant background. In order to avoid a possible 
contribution from — > (pir^, the fit was not extended below 1.895 GeV. The fit yielded 
339 ± 48 D s mesons. The mass value obtained was Mo a = 1967 ± 2 MeV, in agreement 



with the PDG value [18]. The width of the signal was od s = 12.5±1.9 MeV, in agreement 



with the MC estimation. 

A clear signal is seen in the M(K + K-) distribution (Fig. 0c) for the D s region, 1.94 < 
M(K + K^tt ± ) < 2.00 GeV. The fit function for the was a relativistic P-wave Breit- 



Wigner with variable mass and a fixed full-width of 4.43 MeV |L8|], convoluted with a 
Gaussian function whose width, a^, was a free parameter of the fit. The background 
was parametrised with the functional form a[M(K + K~) — 2rrix} b , where is the 



mass. The fit yielded M(0) = 1019.5 ± 0.3 MeV, in agreement with the PDG value WE 



and = 1.7 ± 0.4 MeV, in agreement with MC estimation. The number of mesons 
originating from D s decays was estimated by a side-band subtraction, and was found to 
be in good agreement with the number of D s obtained from the above M(K + K^tt^) fit. 
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Table 1: The differential cross sections da/dp ± a for \rj Da \ < 1.5 and dajdrf" for 3 < 
p® s < 12 GeV for the photoproduction reaction ep — > D S X in various p^ a and rf" bins 
for Q 2 < 1 GeV 2 and 130 < W < 280 GeV. The pf" points are given at the positions 
of the average values of an exponential fit to the p® s distribution in each bin. Statistical 
and systematic uncertainties are quoted separately; the third error corresponds to the 
uncertainty in the D s — > (jnt branching ratio. 



p® s range [GeV] 


(tf s ) fit [GeV] 


da/dp^ B ± stat. ± syst. ± br. 


[nb/GeV] 


3 4 

4 6 
6 12 


3.46 
4.86 
8.09 


2.49 ±0.56 t° il * 
0.523 ± 0.127 ±jj;{jjjf 
0.079 ± 0.021 t° ° l 7 j 


> ± 25% 


7] Da range 




da/drj Ds ± stat. ± syst. ± br. 


[nb] 


-1.5 - -0.5 
-0.5- 0.5 
0.5- 1.5 




1.03 ±0.38 toil ' 
1.90 ±0.35 toil 
0.92 ±0.29 igli? > 


> ± 25% 



Measurement of Inclusive D s Cross Sections 



The inclusive D s cross section is given by: 



where N Ds is the fitted number of D s mesons, C is the integrated luminosity, B Ds ^^ K + K - )ir 
= 0.0177 ± 0.0044 is the combined D s — > (<p-^K + K~) tt decay branching ratio and A is 
the acceptance calculated with the MC sample (Section 3). The MC sample contains 
all events with the D s — > (<f)^K + K~)ir decay channel as well as small admixtures from 
other D s decay modes and from other charm particle decays. Thus these contributions 
were taken into account in the acceptance correction procedure. 

The total D s cross section in the kinematic region Q 2 < 1 GeV 2 , 130 < W < 280 GeV, 
3 < p® s < 12 GeV and —1.5 < r] Ds < 1.5 was measured to be <J ep ^D a x = 3.79 ± 
0.59 (stat.) t oil (syst-) ± 0.94 (br.) nb, where the last error is due to the 25% uncertainty 
in the branching ratio Bp^^. The differential cross sections da/dp® s and da/drj Ds 
are given in Table 1. In Fig. ^| they are compared with the distributions for the D* in 
the same kinematic region. The overall normalisation uncertainty due to the luminosity 
measurement (1.7%) is not included in the cross section errors. 



4.1 Systematic Uncertainties 

A detailed study of possible sources of systematic uncertainties was carried out for all 
measured cross sections by shifting the nominal analysis parameters as described below, 
taking into account resolution effects. For each variation, except for the first one, which 
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is due to the fit systematics, the D s mass and width were fixed to the values found for 
the nominal cut values. The following sources of systematic errors were considered: 



• the uncertainties in the determination of the number of D s mesons were estimated by 
using a quadratic polynomial function for the background parametrisation instead 
of an exponential one and by varying the range of the K + K~tt ± mass distribution 
in the fit procedure; 

• to estimate the uncertainties in the tracking procedure, the track selection cuts, 
including M(K + K~) and cos 6* K , were shifted by at least the expected resolutions 
from the nominal values (Section 3); 

• the cut on p ± s /£J^ >10 ° was changed by ±10%; 

• the MC simulation was found to reproduce the absolute energy scale of the CAL to 
within ±3% [pCfl . The corresponding uncertainty in the cross section was calculated, 
including a shift of 3% in the FLT CAL energy thresholds (Section 3); 

• the fraction of resolved photon processes in the PYTHIA MC sample was varied 
between 40% and 60%; 

• the dE/dx likelihood cuts were changed in the range 0.02 < Z« < 0.10 and 0.08 < 
Li < 0.15; 

• the uncertainty associated with the correction to the true W range was determined 
by moving the Wjb boundary values by the estimated resolution of ±7%. 



None of the above was dominant in the total systematic uncertainty to the inclusive D s 
cross section. All systematic errors were added in quadrature, yielding a total uncertainty 
of ii2%) compared to a statistical error of 16%. 



5 Comparison with QCD Calculations 

The D s cross sections were compared to two types of pQCD calculations || ||. The frac- 
tions of c quarks hadronising as D* and D s mesons were used as input to each calculation. 
The values f(c -> D*+) = 0.235±0.007 (± 0.007) and f(c -> D+) = 0.101±0.009 (± 0.025) 
were extracted pi] by a least-squares procedure from all relevant existing e + e~ experimen- 



tal data f22fl . The errors in brackets are due to uncertainties in the charm hadron decay 
branching ratios. They affect experimental and theoretical cross-section calculations in 
the same way and can be ignored in the comparison. 

The NLO calculation of charm photoproduction in the fixed-order approach of Frixione 
et al. assumes that gluons and light quarks (u, d, s) are the only active partons in 
the structure functions of the proton and the photon. In this approach there is no ex- 
plicit charm excitation component, which can be important in charm photoproduction 
at HERA pj, and charm is only produced dynamically in hard pQCD processes. This 
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calculation is expected to be valid when the c-quark transverse momentum, p±, is not 
much larger than the c-quark mass, m c . 



The structure function parametrisations used in the NLO calculation were MRSG [JHJ for 
the proton and GRV-G HO |L7| for the photon. The renormalisation scale used was fiR = 
m± = \/ml + p]_ } and the factorisation scales of the photon and proton structure functions 
were set to fip = 2m_|_. The pole mass definition is used in this calculation for the c-quark 
mass with a nominal value m c = 1.5 GeV. The Peterson fragmentation function |23|] was 
used for charm fragmentation in this calculation. The Peterson parameter e = 0.035 was 
obtained for D* in a NLO fit @ to ARGUS data |2|. A recent NLO fit ||] to ARGUS 



data yields an e(D s ) value equal to e(D*) within the fit uncertainties. Using the same value 
for both channels leads to the same predictions, except for the difference in f(c — > D* + ) 
and /(c — > which enter the calculation as scale factors. The NLO prediction for the 
total inclusive D s cross section of 2.18 nb is smaller by ~ 2 standard deviations compared 
to the measured cross section. Scaling the Peterson parameter e with the squared ratio 
of the constituent quark masses |23|, m s = 0.5 GeV and m u ^ d = 0.32 GeV [p6fl , leads to 



e(D s ) = 0.085, which yields a NLO prediction 22% lower than that with e(D s ) = 0.035. 

In Fig. ^, the NLO calculations are compared to the differential cross sections. The thick 
curves correspond to the nominal values of and m c , as defined above. For the thin 
curves, a rather extreme value for the c-quark mass, m c = 1.2 GeV, and a fiR value of 
0.5m± have been used. The D s cross section decreases steeply with rising p^, as in the 
D* case. The NLO calculation reproduces within errors the shape of the p^ s distribution 
but underestimates the data for the nominal parameter set. For the r] Ds distribution the 
NLO predictions are below the data in the central and forward regions. A similar effect 
was observed when da/dr] D distributions were compared with various NLO predictions 
over a wide range of W values and photon virtualities ^ Q . 



Recently Berezhnoy, Kiselev and Likhoded (BKL) have suggested another model which 
does not employ any specific fragmentation function. In this tree-level pQCD 0{aa^ s ) 
calculation, the (c, q) state produced in pQCD is hadronised, taking into account both 
colour-singlet and colour-octet contributions. Using the experimental value for f(c — > 
D* + ), a c-quark mass m c = 1.5 GeV, a light constituent quark mass m q = 0.3 GeV and 
the CTEQ4M proton structure function parametrisation |[29|| , the ratio of the colour-octet 
and colour-singlet components was tuned to describe the ZEUS D* photoproduction cross 
sections This ratio and the experimental value for f(c — > Df) were used to obtain 
predictions for D s photoproduction. In this strange quark mass, m s = 0.5 GeV, 

was used instead of m q . The calculated cross sections for D s and -D s * ± were combined, 
since the inclusive D s channel includes fully the prompt -D s * ± meson production. 

In Fig. |3|, the BKL calculations are compared to the D s differential cross-section mea- 
surements. The agreement with the data is better than that of the NLO calculation 
with the nominal parameters, but the shape of the 7] Ds distribution is not matched by 
the BKL prediction. The total predicted D s cross section in the kinematic range of the 
measurement (3.37 nb) is compatible with the measured inclusive cross section. 



5 A comparison of the D* data with the BKL calculations can be found in the BKL paper ||. 
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6 D s to D* cross-section ratio and 7 S 

In the D s kinematic region, as defined in Section 4, the D* cross section was measured to 
be a ep ^D*x = 9.17±0.35 (stat.) 1^39 ( s Y s t-) n b 0- This yields a ratio (y ep -*D s x I '(Jep^>D*x = 
0.41±0.07 (stat.)lo!o5 ( s Y s t-) ±0.10 (br.), where common systematic errors (Wjb and CAL 
energy scale) have been removed and the last error is the uncertainty in Bd^^. This 
ratio is in good agreement with the ratio f(c — > Df)/ f(c — > D* + ) = 0.43±0.04±0.11 (br.) 
obtained from results of e + e~ experiments (see Section 5). The last error originates from 
the uncertainty in Bo^^n and can be ignored in the comparison. 

The strangeness-suppression factor, j s , is the ratio of probabilities to create s and u, d 
quarks during the fragmentation process. In simulation programs based on the Lund 
string fragmentation scheme PB|, 7 S is a free parameter with a default value of 0.3. By 



varying the value of 7 S in the PYTHIA simulation [|TJ|] , a direct relation can be obtained 



between 7 S and the D s to D* cross-section ratio. Fixing the value of f(c — > D* + ) in 
PYTHIA to 0.235 yields 

7, = 0.27 ± 0.04 (stat.) t^s ( s Y st -) ± °- 01 ( frac ± °- 07 ( br • 

The third error is due to the uncertainty in f(c — > D* + ), while the forth one results from 
the uncertainty in B Ds ^ n . Adding all errors in quadrature, except the last one, gives 
7 S = 0.27 ± 0.05 ± 0.07 (br.). 

Previously, 7 S was measured mainly from the ratio of K to 7r production and from the mo- 



mentum spectrum of K mesons in hadron-hadron [|30], |3l| and e + e~ [[32], |33|, [34| collisions, 
as well as in deep inelastic scattering (DIS) experiments |35], |3B|, |57| . The most accurate 
measurement, obtained in pp collisions |3l| , is 7 S = 0.29 ± 0.02 (stat.) ± 0.01 (syst.). The 
DIS results require a lower value, 7 S ~ 0.2. Recent results from e + e~ collisions are in some 
disagreement with each other. The SLD preliminary result [B3j is 7 S = 0.26 ± 0.06, while 



OPAL finds 7 S = 0.422 ± 0.049 ± 0.059. Previous 7 S measurements are therefore in 
good agreement with that of this analysis, except the latest OPAL value, which is about 
2 standard deviations higher. 

Existing 7 S values obtained from heavy-meson production (charm and beauty) in e + e~ 
collisions |32] are centred around 0.3. For charm production, the ratio 2/(c — > Df)/[f(c — * 
D + ) + f(c — > D )] was used as a measure of 7 S . Using for the above fractions the 
more recent values quoted in pl| leads to 7 S = 0.26 ± 0.03 ± 0.07 (br.), where the latter 



uncertainty originates from Bp^^ and can be ignored in the comparison with the ZEUS 
result. 

The results presented here on the D s to D* cross-section ratio and on 7 S , taken together 
with charm production data in e + e~ annihilation, tend to support the universality of 
charm fragmentation. 



7 Summary and Conclusions 

The first measurement at HERA of inclusive photoproduction has been performed 
with the ZEUS detector. The cross section for Q 2 < 1 GeV 2 , 130 < W < 280 GeV, 
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3 < p1 a < 12GeV and -1.5 < rj D ' < 1.5 is cr ep ^ DsX = 3.79 ± 0.59 (stat.) 1 ^ ( s y st ± 
0.94(br.)nb. The differential cross sections da /dp® 3 and da/drj Ds are generally above 
fixed-order NLO calculations, as was the case with the results previously obtained for D* 
photoproduction in the same kinematic region. The BKL calculation, using the octet- 
to-singlet ratio tuned to fit the ZEUS D* data, describes the D s cross sections reason- 
ably well, but the shape of the r] Ds distribution is not matched by the BKL prediction. 
The cross-section ratio, a ep ->D s x/o'ep-*D*x, in the kinematic region as defined above is 
0.41 ± 0.07 (stat.)^;^ (syst.) ± 0.10 (br.), in good agreement with the ratio of c quarks 
hadronising into D s and D* mesons, extracted from e + e~ experiments. From this ratio, 
the strangeness-suppression factor in charm photoproduction, within the LUND string 
fragmentation model, has been calculated to be 7 S = 0.27 ± 0.04 (stat.) (syst.) ± 
0.01 (frac.) ± 0.07 (br.), in good agreement with the 7 S value extracted from charm pro- 
duction in e + e~ annihilation. 
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Figure 1: (a) Distribution of the sample used for fitting the dE/dx parametrisation in 
the (dE/dx, p) plane. Here dE/dx is normalised to a minimum ionizing particle (MIP), 
defined as the average truncated mean of pion tracks in the momentum range 0.3 < p < 
0.4 GeV. The lines indicate the fit result [HJ. (b) M(K + K~tt ± ) distribution for events 
inside the <fi mass range (1.0115 < M(K + K~) < 1.0275 GeV). The solid curve is a fit to 
a Gaussian (representing the resonance) plus an exponential background, (c) M(K + K~) 
distribution for events inside the D s mass range (1.94 < M(K + K~tt ± ) < 2.00 GeV). The 
solid curve is a fit to a Breit-Wigner convoluted with a Gaussian-shaped resonance and a 
background parametrisation of the form a[M(K + K~) — 2m^] fc . 
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Figure 2: Differential cross sections for the photoproduction reaction ep — > D X: 
(a) da/dpf and (b) da/dr/ D , where D stands for D* or D s . The pf points are drawn 
at the position of the average value of an exponential fit in each bin. The r] D points are 
drawn at the middle of each bin. The inner error bars show the statistical uncertainty, 
while the outer ones show the statistical and systematic errors added in quadrature. Nor- 
malisation uncertainties due to the D s — > 0tt branching ratio are not included in the 
systematic errors or in the theoretical calculations. The D s (dots) and D* (triangles) 
data are compared with the NLO predictions for D s (full curves) and D* (dashed curves) 
with two parameter settings: m c = 1.5 GeV, /Ir = m± (thick curves) and m c = 1.2 GeV, 
Hr = 0.5m_|_ (thin curves). 
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Figure 3: Differential cross sections for the photoproduction reaction ep — > D S X: 
(a) da/dpZ and (b) da/drj D " compared to the BKL model ||. Colour-singlet (dashed 
curves) and colour-octet (dotted curves) contributions are plotted separately. Their sum 
is shown as the full curves. The p_^ s points are drawn at the position of the average value 
of an exponential fit in each bin. The r] Ds points are drawn at the middle of each bin. 
The inner error bars show the statistical uncertainty, while the outer ones show the statis- 
tical and systematic errors added in quadrature. Normalisation uncertainties due to the 
D s — > <f)ir branching ratio are not included in the systematic errors or in the theoretical 
calculations. 
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